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A formalism for hadron production at high pr in heavy-ion collisions has been developed such 
that all partons hadronize by recombination. The fragmentation of a hard parton is accounted for by 
the recombination of shower partons that it creates. Such shower partons can also recombine with 
the thermal partons to form particles that dominate over all other possible modes of hadronization 
in the 3 < pr < 8 GeV range. The results for the high pr spectra of pion, kaon, and proton 
agree well with experiments. Energy loss of partons in the dense medium is taken into account 
on the average by an effective parameter by fitting data, and is found to be universal independent 
of the type of particles produced, as it should. Due to the recombination of thermal and shower 
partons, the structure of jets produced in nuclear collisions is different from that in pp collisions. 
The consequence on same-side correlations is discussed. 



I. INTRODUCTION 

In the production of hadrons at high px in heavy-ion 
collisions there are by now three theoretical collabora- 
tions that have shown the importance of quark recombi- 
nation 0, 0, • While there are some differences among 
the three approaches, they all agree in the basics and 
in the successful interpretation of the experimental data, 
among which the most outstanding ones are the p/ir ra- 
tio being around 1 in the 3 < px < 4 GeV range Q 
and the scaling law of elliptic flow in the number of con- 
stituents pj, Lj • The differences concern mostly with the 
treatment of hadronization in the 4 < px < 8 GeV range. 
In this paper we show how the particles produced in that 
range arise from the recombination of the thermal par- 
tons and the shower partons created by hard partons at 
high pt- The determination of the shower parton dis- 
tributions (SPD) has recently been achieved by studying 
the fragmentation functions (FF) in the framework of the 
recombination model 0. This paper contains the first 
application of the SPD outside the realm of parton frag- 
mentation. As a consequence we find a new component 
that stands between the recombination of soft thermal 
partons at lower px and the fragmentation of hard par- 
tons at higher px |3j , and is also different from the direct 
recombination of soft and hard partons • 

In we have avoided the need to specify the origin of 
the partons before recombination in order to be indepen- 
dent of the models describing the early and intermediate 
phases of the evolution in heavy-ion collisions. We use 
the measured pion spectra to infer the quark and anti- 
quark distributions just before hadronization, and then 
use those distributions to determine the proton spectrum. 
In that way the calculated p/it ratio is a direct conse- 
quence of the recombination model with essentially no 
dependence on the other aspects of the dynamics such 
as the separation into soft and hard components. In 
this paper we do enter into the origins of the partons. 
The main difference between our approach and the con- 



ventional treatment of particle production at high px is 
that we do not use FF to represent the hadronization 
of hard partons produced at higher px- In our view all 
hadrons are produced by recombination at any px- FF's 
are phenomenological functions that do not specify the 
hadronization mechanism. Although string models can 
be useful in the description of hadronization of a pair of 
q and q receding from each other after being created in 
vacuum in, for example, e + e~ annihilation, they cannot 
be applied to heavy-ion collisions where the abundance 
of color charges renders invalid any notion of stretched 
color flux tubes between pairs of partons produced in 
hard collisions . In the absence of such models to give 
a conceptual basis for fragmentation, it is necessary to 
find a meaningful hadronization scheme for hard partons 
at high px- 

Our view is that a hard parton creates a shower of 
partons that recombine subsequently to form hadrons. 
Although the branching process of gluon radiation and 
pair creation that eventually lead to shower partons at 
low virtuality cannot be calculated, the distributions 
of the shower partons can nevertheless be determined 
from known FF's in the framework of the recombination 
model, analogous to how the q and q distributions are de- 
termined from the experimental distribution dN/pxdpx 
of pions, as done in 0. In @ a variety of SPD's are given 
as functions of the momentum fractions of the shower 
partons. In principle, there should be dependences on 
Q 2 ; however, for use in heavy-ion collisions at RHIC, the 
SPD's given in for Q 2 = 100 GeV 2 is adequate, since 
the dependence on Q 2 is not severe. 

With the SPD's at hand, we can now consider their 
role in heavy-ion collisions. Hard scattering of partons 
gives rise to partons with high px, which undergo energy 
degradation as they traverse the dense medium Il0| . 
Instead of tracing the spatial coordinates of the hard par- 
tons, we shall use an effective parameter £ to describe 
the average fraction of partons that escape the dense 
medium and are able to hadronize outside. The value of 
£ will be determined phenomenologically, and will repre- 
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sent an independent check on the degree of jet quenching 
[Tl|. A more important issue is the px dependences of 
the hadrons detected and the origins of the partons that 
contribute to the formation of those hadrons at different 
regions of px ■ The recombination of thermal and shower 
partons will be shown to be important with the conse- 
quence that the structure of jets produced in heavy-ion 
collisions is different from that produced in pp collisions. 
That difference will manifest itself in the correlation be- 
tween particles in the same jet. 

The recombination of two shower partons in the same 
jet is the same as the usual fragmentation of hard partons 
and becomes important at higher px- The recombination 
of two shower partons arising from two neighboring hard 
partons is also possible, but will not become important 
until the collision energy is very high, e.g., at LHC. 

In what follows the production of mesons (pions and 
kaons) will be considered first, and then the baryons 
(specifically proton). Same-side correlation will be dis- 
cussed, but only qualitatively to explain some apparent 
puzzle in the data. 



II. MESON PRODUCTION BY 
RECOMBINATION 

A. General considerations 

In an overview of a heavy-ion collision we can divide 
the process into two stages: (I) the initial evolution- 
ary phase, and (2) the final hadronization phase. Our 
concern in this paper is mainly on the latter phase, al- 
though what partons recombine depends on the former. 
A proper formulation of the formation of a dense medium 
must be done in full 3D-space and time, but to de- 
scribe hadronization at large transverse momentum can 
be much simpler. Since the spatial volume in which re- 
combination can occur is small, any two partons that 
are not collinear cannot recombine. Thus collinear par- 
ton momenta in ID are all that we need to consider for 
the production of a particle with a particular momentum 
p. The spatial extension along that direction is also con- 
strained to the size of the hadron, so an integration of the 
remaining spatial variable results in a momentum-space 
description in which the parton momenta should reflect 
the wave function in momentum-space representation of 
the produced hadron. We can therefore formulate the 
recombination process in the ID momentum space only. 
Starting with p we investigate the partons that can re- 
combine to form the particle at that p. Of course, by 
not starting with a 3D space from the outset, we can- 
not calculate the density and expansion properties of the 
soft thermal partons, which will have to be introduced 
phenomenologically. A more elaborate formulation can 
be given in terms of Wigner functions in 6D coordinate- 
momentum space 0, Q ■ Such formulations do provide a 
more complete picture of the collision process with the 
added advantage of being able to specify the soft compo- 



nent dynamically. It should be emphasized that whereas 
the recombination process is considered by us in ID, it 
does not mean that the formalism is insensitive to the 
realistic problem of heavy-ion collisions and the 3D na- 
ture of the colliding nuclei. We shall return to this point 
when we discuss centrality dependence. As far as the 
hadronization part of the problem is concerned, our for- 
mulation in the ID momentum space is consistent with 
the result of the formulation in 6D coordinate-momentum 
space after integration over the other 5 variables. 

In our present problem of heavy-ion collisions we shall 
consider particle production in the transverse plane at 
rapidity y — only. Although the transverse plane is 2- 
dimensional, we shall consider only the direction in which 
a hadron is detected so that all partons relevant to the 
formation of such a particle move in the same direction. 
Thus it is unnecessary to carry the subscript T on all our 
momentum labels to denote "transverse." The invariant 
phase space element is therefore dp/p°, which we shall 
approximate by dp/p for relativistic particles. 

In the recombination model [l2T ] the invariant inclusive 
distribution for a produced meson with momentum p is 

dN M f dp 1 dp 2 

V—, — = / F qi j,(p 1 ,p2)RM{pi,P2,P) , (1) 

dp J Pi P2 

where F q qi(pi,p2) is the joint distribution of a quark q 
at p\ and an antiquark q' at p2, and Rm(j>i,P2,p) is the 
recombination function (RF) for qq 1 — > M. For central 
collisions it is immaterial which direction p points. Al- 
though we only consider p in the transverse plane here, 
Eq. Q has been used in the longitudinal direction in 
0[b|. where the RF is specified. It is 

I / \ ( \ 

= WTYJTT) \i) (?) 

.*(* + *-lY (2) 
VP P J 

where B(m, n) is the beta function. For pion it is shown 
in 13] from the analysis of Drell-Yan production data in 
pion-initiated process |14j in the framework of the valon 
model 12] that a = b = 0. Thus we have explicitly 

Rn(Pl,P2,P) = — o-<M 1 1 , (3 

P 2 \P P J 

The statistical factor for the recombination process is 1 
0] . For kaon it follows from the constituent quark masses 
of the valons and kaon-initiated inclusive production [HH ] 
that a = 1 and b — 2. 

The (^-function in Eq. @ guarantees the conservation 
of momentum in the recombination process. Rm is an 
invariant distribution that is related to the non-invariant 
probability density Gm (j/i : 2/2 ) of finding the two valons 
in M with momentum fractions y\ and y 2 by 

Rm(pi,P2,p) = y\y2G M {yiy2), y% = Pi/p ■ (4) 
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Thus for pion (2/13/2) is a constant, apart from the 5- 
function. It means that in the momentum space the wave 
function of pion in terms of the valons is very broad, 
corresponding to the pion being a tightly bound state of 
its constituent quarks. That is not the case for kaon. The 
broadness of G 7r (xiX2) is an important reason why the 
thermal-shower recombination in the formation of pions 
makes a dominant contribution in the intermediate px 
range, as we shall see below. 

It should be noted that given in Eq. © differs 
from the RF given in , where our formulation of the 
high-pT problem is in the 2D transverse plane. Even 
the dimension of in [l| is equivalently pTp 2 , whereas 
it is dimensionlcss in Eq. Thus caution should be 
exercised in relating our treatment here to that in 

It is important to recognize that the RF describes the 
probability of recombination of quarks (and antiquarks). 
In the recombination model the gluons are not regarded 
as partons that can directly hadronize [l2l Il3| . They 
hadronize through the q and q channels by conversion 
to qq pairs. Thus F qq i in Eq. Q must include all the q 
and q' generated by gluon conversion for the purpose of 
hadronization through the use of Rm in Eq. . The sea 
is therefore saturated by complete conversion from the 
gluons, a procedure that leads to the correct inclusive 
cross section for the hadronic production of pions, both in 
normalization and in the momentum distribution Il2l 
[l3]| . The saturation of the sea of q and q' makes possible 
that the momenta of the gluons (which carry roughly 
half of the momentum of each nucleon) can be properly 
accounted for in the produced hadrons in any hadronic 
or nuclear collision. This procedure of saturating the sea 
is also followed implicitly in for the determination of 
the shower partons in the fragmentation of an initiating 
parton. We also note that the conversion of gluons to 
quark pairs increases the number of degrees of freedom of 
all the partons in the medium and thereby overcomes the 
decrease of entropy that one might otherwise conclude by 
focussing on individual recombination processes. 

Lastly, we mention that soft gluon emission and ab- 
sorption are always possible in a recombination process; 
in fact, that is how color mutation can take place for a qq' 
pair to become colorless, while the partons dress them- 
selves to become valons before forming a meson. Such 
soft processes do not change the momenta of valons and 
therefore do not influence the momentum consideration 
in Eq. @ ■ We do not consider quark-antiquark-gluon re- 
combination with large momentum fraction for the gluon 
because the valon representation is complete. That is 
precisely the reason why the valon model was constructed 
in the beginning to describe hadron structure on the one 
hand as well as the recombination function (for the time- 
reversed process) on the other 0] ■ 



B. Pion distribution 

Restricting our attention to pion production for now, 
we obtain from Eqs. and @ 

dN n IF, 

— — = -3 / dpiF^ijpx^- pi), (5) 
pdp p 6 J 

which clearly exhibits the simple dependence of the pion 
spectrum on the momentum distributions of the q and q' 
that recombine. It should be recognized that, since we 
work in ID, dN^/pdp in Eq. (J5J is actually the number 
density of pions in pTdp^dydfy evaluated at y = 0. With 
the assumption that it is independent of 4> in central col- 
lisions, what we denote as dN^/pdp is equivalent to the 
experimental dN / (2TrpxdpT) , where the number N refers 
to the integrated result over all </> so that when divided 
by 2ir the average density in 2D is the same as what we 
calculate in ID. 

Note that in Eq. J5J Pi is integrated over a wide range 
due to the broadness of the RF. There are two compo- 
nents of parton sources that contribute to F qq -/ . One is 
thermal (T) and the other shower (S). Before giving 
the specifics of what they are, let us first express F qq > in 
terms of them in a schematic way: 

F q q' =TT + TS + (SS), + (SS) 2 . (6) 

All four terms make contributions at all pr, although 
each is important in only restricted regions of pr- {SS)i 
denotes two shower partons arising from one hard parton 
(hence within one jet), and can be related through the 
use of R to the usual fragmentation that is described by 
the D function. (SS)?, denotes two shower partons that 
are from two separate but nearby hard partons, and are 
therefore associated with two overlapping jets. (SS)2 is 
not expected to be important unless the density of hard 
partons is extremely high, such as that possibly at LHC. 
For brevity, we shall write (SS) 1 simply as SS, when 
no confusion is likely to arise. TT signifies two thermal 
partons whose recombination yields the thermal hadrons, 
usually referred to as the soft component. TS denotes 
thermal-shower pairing and is the new component that 
has never been considered before. It turns out to be 
important in the 3 < pt < 8 GeV range. We emphasize 
that the TS term would be absent if we do not treat 
the fragmentation of a hard parton as the recombination 
of shower partons as done in Q • It is now evident that 
by considering S as the showering effect of hard partons 
all hadrons are produced by recombination, as is made 
explicit by Eqs. JSJ and (j^Jl in the case of pions. Let us 
now specify what T and S represent. 

1. Thermal component 

T is the thermal component including hydrodynamical 
flow. It is not our intention to derive soft parton distribu- 
tion from hydrodynamics. We shall simply assume what 
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is necessary to give rise to the observed distribution of 
pions for p T < 2 GeV. Since the observed dN^/pdp at 
low p is exponential, the necessary invariant parton dis- 
tribution for the thermal component is 



dN th 

T[pi) = pi— — = Cpi exp(-pi/T), 
dpi 



(7) 



where T is the inverse slope enhanced by flow. C is the 
normalization factor to be adjusted to fit the pion data 
at low p, and has the dimension [momentum]" 1 . For 
noncentral collisions, which we do not consider in this 
paper, C would depend on centrality. We assume that 
the thermal partons are uncorrelated, except in special 
circumstances, so that we may write F^, in the factoriz- 
able form 

F$(pi,P2) = T(pi)T( P2 ) = C 2 Pl p 2 eM~{Pi +P2)/T] .(8) 

Substituting this into Eq. JSJ yields for the thermal com- 
ponent of the pion distribution 



regarded as an empirical quantification of the degree of 
energy loss. 

The hard parton i at momentum k creates a parton 
shower. Since there are various types of shower partons 
for each type of initiating parton i, let us use Sj to denote 
the matrix of SPD's for i — > j. Although i can be u, d, 
s, u, d, s and <?, j is allowed to be quark and antiquarks, 
but not gluon. The role of gluons in the recombination 
process has already been discussed earlier at the end of 
Sec. II-A. The sea partons in the shower are saturated 
by gluon conversion. In the notation of we use Kns 
to denote the valence quark in the shower, L (L s ) the 
light (strange) sea quarks in a quark-initiated shower, 
and G (G s ) the light (strange) sea quarks in a gluon- 
initiated shower. Thus the shower matrix Sf has the 
form 



dK h G 2 , 
^^expHVT), 



which is the exponential form aimed for. To fit the data 
the parameters are 



C = 23.2 GeV" 1 , 
as we shall see below. 



T = 0.317 GeV, (10) 



2. Shower partons 

Next we consider the shower distribution. In a heavy- 
ion collision let the parton i be scattered into the trans- 
verse plane at rapidity y = with probability fi(k), i.e., 



dN, 



hard 



kdkdy 



fi{k) 



(11) 



y=o 



where k is the transverse momenta of the parton. We 
shall use the parametrization of fi(k) given in Ref. |15| . 
obtained for the study of dilepton production in central 
collision of gold nuclei at ^/snn — 200 GeV. Due to en- 
ergy loss of the partons in the dense medium, not all 
partons emerge from the reaction zone to hadronize out- 
side. Only a fraction of them do, and we shall use £ 
to represent the effective fraction after averaging over 
all central events such that denotes the number 

of unquenched partons with momentum k that are to 
hadronize. Note that we do not refer to them as jets, 
since the notion of jets presupposes that a hard parton 
fragments into a jet of hadrons. That supposition is, of 
course, just what we want to avoid. We regard £ as an 
effective fraction because we do not consider its depen- 
dence on k and do not delve into the space-time proper- 
ties of the hard scattering and subsequent evolution. The 
value of £ will be determined phenomenologically, and be 



s- 
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L s 
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K s 


Kg 


G 


Gs 



= u,d,s,g, j = u,d,s, 



(9) where K = K NS 



(12) 

L and K s = K^s + L s . The anti- 
quarks u, d and s have the same structure, and are related 
to u, d, s as sea, and vice-versa. The parametrizations 
for these SPD's have been completely determined in 
as functions of the momentum fraction z of parton j in 
parton i. Note that in Eq. (|12fl there is no fourth column 
corresponding to gluons in the shower. The determina- 
tion of the five essential SPD's is carried out on the con- 
dition of no shower gluons so that all hadrons specified 
by the fragmentation functions are formed by qq' recom- 
bination without gluons. The underlying physics for this 
has already been discussed in the last two paragraphs of 
Sec. II.A. 

The distribution of shower parton j with transverse 
momentum p\ in central heavy-ion collisions is then 



r-OQ 

S(Pi) = tJ2 dkkMk)Si( Pl /k) 

J kn 



(13) 



where p\ and k are collinear. Since the input on hard 
parton distribution /,(fc) cannot be valid at low k, we 
shall consider the above integral only for k > fcoj f° r 
which we set the minimum at ko = 3 GeV. In practice 
we cut off the upper limit of integration at 20 GeV. 



3. Thermal-shower recombination 

With the shower partons specified we can now combine 
them with the thermal partons to describe the TS term 
in Eq. ©. We have 

T{p 1 )S( P2 )=iCp 1 e-^l T J2 J dkkf t (k)Si(p 2 /k) , 

where the distributions of thermal light quarks are as- 
sumed to be flavor independent and the appropriate one 



5 



is implied to pair off with j to form the meson under con- 
sideration. The contribution to the pion spectrum from 
thermal-shower recombination is then, using Eq. 



pdp 



1 f p 

— / dp 1 T(p 1 )S(p-p 1 ) 
P A Jo 



(15) 



where a sum over j is implied to match the flavors of 
the valence quarks of the detected pion. Only the overall 
normalization of this term depends on £. The p depen- 
dence is our prediction, which is used to fit the data and 
thereby determine £. 



J 



4- Shower-shower recombination 



For two shower partons in the same jet we have 



1>2 



Pi 



where the curly brackets signify the symmetrization of 
the leading parton momentum fraction 



GL6) 



•~2 



1-21 



sHzi)si 



Z-2 



1-Zi 



Si 



1 



Z-2 



S 3 Az 2 ) 



(17) 



We have shown in that the SPD's can be determined from the recombination formula for the fragmentation function 

•i'2 



dx\ dx2 

X\ X 2 



1 — Xi 



R M {xi,x 2 ,x) . 



(18) 



r 



A sum over j and j' is implied in consort with the j and 
j' labels hidden in the RF that are relevant for M. The 
substitution of Eq. (|16fl for the (SS) i term in Eq. © 
into Eq. JQ) clearly yields 



(19) 



which is the usual formula for the production of a meson 
at high pt in the fragmentation model, except for the 



J 



presence of £ here for reasons that have already been 
discussed above. 



There is finally the possibility of recombination of two 
shower partons from two different but partially overlap- 
ping showers. The corresponding {SS)2 term in Eq. © 
should then be 



(SS) 2 ( Pl ,p 2 ) = 6 V ^ 2 Y, J dkdk'kk'fi{k)f v {k'\ 



(20) 



where a multiplicative factor 5 y $ is included to reflect 
the probability of overlap in y and 4> of the two showers 
in order for collinear recombination of the partons j and 
j' to take place. The value of 5 y( f, can be estimated by 
studying the size of the jet cone, and is expected to be 
small. Thus this mode of recombination is not likely to be 
important at RHIC. However, at very high energy, such 
as at LHC, where f%{k) is orders of magnitude higher, 
the (SS) 2 term may well become significant. 



5. Result on the pion spectrum 

Collecting all the pieces of Eq. © together and substi- 
tuting them in Eq. we obtain the four contributions 
to the pion spectrum. The parameters C and T in Eq. 
(|10f) are determined by fitting the low-px data. Ignor- 
ing the (SS) 2 contribution on the basis that 5 y( j, is very 
small, there is only one parameter, £, to adjust to fit the 
data for px > 2 GeV. The result is shown in Fig. 1. With 
the value 



£ = 0.07 



(21) 
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FIG. 1: Transverse momentum distribution of it in Au-Au 
collisions. Data are from |lo| . The solid line is the sum of four 
contributions to the recombination of partons: TT (dashed 
line); TS (line with crosses); (SS)i, two shower partons in 
one jet (line with open circles); [SS)2, two shower partons 
from two overlapping jets (line with squares). 



the fit of the tt° data from PHENIX [16| on central Au-Au 
collisions at ^/smn = 200 GeV is excellent up to pt — 10 
GeV. Note that in the region 3 < pr < 8 GeV, the domi- 
nant contribution is from thermal-shower recombination 
(line with crosses). The conventional jet fragmentation is 
from shower-shower recombination in one jet (line with 
circles); it becomes more important than the TS con- 
tribution only for px > 9 GeV. To show the relative 
size of the (SS)2 recombination from 2 jets, we assume 
Sytjj — 0.01 (a very rough estimate) just to put its con- 
tribution on the figure. It is indicated by the line with 
squares, which is much lower than all others. 

The result that thermal-shower dominates over shower- 
shower (1-jet) recombination for pr < 8 GeV is our main 
finding in this work. It shows the importance of consider- 
ing the interaction between the thermal partons and the 
partons created by hard scattering. That interaction be- 
comes particularly significant at the hadronization scale 
where recombination occurs. 

One way to see why TS is greater than SS in their 
contributions to Eq. JQ) via © is to examine their sum 
in the following form 



TS + SS = £ 



dkkfi(k)Sf 



Pi 



P2 



(22) 



r 



putting aside the other symmetrizing term in SS with- 
out any impediment to our argument. The first (thermal) 
term inside the square bracket is much larger than the 
second (shower) term when p\ is small (but not infinites- 
imal), which is a region of p\ that is relevant for the soft 
parton to recombine with a shower parton at pi > 3 GeV 
only if R 7r {pi,p2,p) is broad enough to encompass both. 
This linking between soft and semi-hard partons not only 
enhances the spectrum over simple fragmentation, but 
also modifies the structure of what is usually regarded 
as minijet. Since there are no thermal soft partons in pp 
collisions, the pion distribution at intermediate and high 
Pt in pp collisions must differ from that in A A collisions. 
Furthermore, the same-side correlations in the two cases 
are necessarily different, as we shall discuss in Sec. 4. 



6. Dependence on the recombination function 

To see how the thermal-shower recombination depends 
on the broadness of RF, consider the general form of the 
valon distribution 



i 



B(a + l,a + 1 



(yiy 2 ) a %i + y 2 -i).(23) 



Equations © and (0} follow from Eq. (|23J| for a = 0. 

If a were larger, as one would expect for the p meson 



(since small a corresponds to a tightly bound state), then 
the valon distribution would be more sharply peaked at 
Vi =2/2 = 1/2. In that case the recombination of a ther- 
mal parton at low p\ and a shower parton at interme- 
diate P2 is suppressed. To make this point transparent, 
we show in Fig. 2(a) several possible widths of a single- 
valon distribution G(y) obtained from G(yi, j/2) by one 
integration, i.e., 



G(y) = 



1 



B(a +1,0+1" 



(24) 



which is normalized to 1 by one more integration. The 
corresponding RF is given by Eqs. © and @. When 
we use that RF in Eq. and calculate the distribution 
dN/pdp with only the TS contribution to F q q taken into 
account, the result is shown in Fig. 2(b) for a = 0, 1, 2, 5. 
It is evident that the recombination of thermal-shower 
partons is significantly suppressed at high px when a is 
large, and becomes negligible when R(p\,p2,p) tends to- 
ward being proportional to S(pi — p/2)5{p2 — p/2). 



C. Energy loss 

In the absence of a space-time study of the problem 
that includes the locations where hard collisions occur, 
it is not possible to consider the medium effect on each 
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FIG. 2: (a) Various valon distributions in momentum fraction 
y according to Eq. 1241 : 

(b) The corresponding pr distributions from TS recombina- 
tion. 



and every hard parton that traverses the medium. We 
have used the parameter £ to represent the overall effect 
after averaging over all events in central collisions. The 
value £ = 0.07 given in Eq. lj2"T|l is a quantity deduced 
from fitting the pion spectrum. It clearly indicates that 
not all hard partons created in a heavy-ion collisions can 
get out of the dense medium to hadronize. As suggested 
by the STAR data 0, only those near the surface can 
escape the quenching effect. Our result on the value of 
£ seems low by comparison to the nuclear modification 
factor Raa(pt), which is roughly 0.2 for pt > 6 GeV in 
central Au-Au collisions. To understand their difference, 
let us examine what they measure respectively. 
Raa{pt) is defined by the ratio 



Raa(j>i 



dN/p T dp T (AA) 
N coU dN/p T dpT (pp) 



(25) 



where N co u is the average number of binary collisions. 
In the denominator the production of particles at high 
Pt in pp collisions can be well described by the fragmen- 
tation of hard partons. However, we have seen that the 
usual fragmentation corresponds to the shower-shower re- 
combination in 1-jet, which is not the important part of 
hadronization in AA collisions for 3 < px < 8 GeV. The 
suppression factor £ may, in the spirit of Eq. (|25p. be 
written in the schematic form 



l dN/p T dp T {AA) \ 
\ JTSR(AA) / 



(26) 



where S is the shower component expressed in Eq. (|13|1 . 
but without the £ factor. The angular brackets denote 
an average over all px ■ The denominator is the thermal- 
shower recombination in AA collisions, if all hard partons 
get out of the medium to hadronize. It is then clear why 
£ is smaller than Raa- It should not only account for the 



suppressed number of hard partons in the numerator that 
get out from the dense medium to hadronize (which Raa 
does also), but is also made smaller by the denominator 
where the TS recombination is larger than the scaled 
contribution from pp collisions. 




6 P T (GeV) 8 
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FIG. 3: Nuclear modification factor Raa compared to the 
calculated result using only pions, R\a- 

Another way to exhibit the effect of energy loss on 
our result is to calculate Raa{pt) directly from our pion 
distribution, assuming that the pions dominate over all 
other particles. That assumption is invalid for pr < 4 
GeV, since the production of proton is known to be 
roughly equal to that of pion in the 3 < pt < 4 GeV 
region. Nevertheless, it is illuminating to see what the 
pion component gives (call it R A a) relative to the exper- 
imental Raa for pt > 4 GeV. Since SS recombination 
is just the fragmentation component in AA collisions as 
stated in Eq. ifT^jl . we can obtain Nco^dN^ / prdpripp) 
simply by omitting the factor £ from dN^ Jpdp. That 
is, we have 



Raa(p) 



dN^/pdp 
Z^dNSS/pdp 



(27) 



where the numerator corresponds to the solid line in Fig. 
1 and dN^ s jpdp corresponds to the line with open circles 
in the same figure. The result for R\ A ^ s shown in Fig. 3 
by the solid line, which is in reasonable agreement with 
the data on Raa [13 that is not limited to the pions. It 
is now clear that the ratio of those two lines is equal to 
the ratio R\a/^ which in turn is approximately equal to 
Raa/S,- Since the thermal-thermal recombination is neg- 
ligible for pt > 4 GeV, the ratio (dN^/pdp)/(dN^ s jpdp) 
is essentially independent of £. Thus R\aI^ x& roughly 
independent of energy loss and therefore of centrality. 

It is worth remarking on how the nuclear-size effect en- 
ters into our one-dimensional treatment of the hadroniza- 
tion process. We emphasize that hadronization occurs 
outside the volume of dense matter and along the direc- 
tion of the detected particle. The 3D nuclear-size effect 
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influences the properties of the partons before they enter 
into the interaction region for hadronization, and is there- 
fore not explicitly present in the recombination formula. 
The recombining partons are either the soft thermal ones 
or the semi-hard shower partons. Centrality obviously af- 
fects the magnitude of the thermal source, i.e., C, and 
the degree of energy loss, i.e., £. We have made prelimi- 
nary study of the centrality dependence of dN^/pdp and 
found agreement with the data by appropriate adjust- 
ment of C and £ in essentially the same way as we have 
done in this paper for the most central collisions. We 
mention this only to point out that the ID description 
of recombination does not preclude the possibility of ac- 
counting for the nuclear-size effect, which is manifestly 
3D, but occurs prior to the hadronization process. 

The suppression factor £, as expressed in Eq. <|26|l . can- 
not be measured directly. There is, however, a measur- 
able quantity that is closely related to £. Since only the 
hard partons that are created near the surface of the 
collision region can get out, most jets detected in AA 
collisions do not have a partner in the opposite direction. 
That does not mean the non-existence of back-to-back 
jets. A pair of hard partons that are created near the 
edge, but directed at around 90° relative to the radial 
position of creation measured from the center, and yet 
remain in the transverse plane (i.e., roughly tangent to 
the cylinder), do not go through the bulk of the medium. 
Those two hard partons can then lead to back-to-back 
jets that are detectable. It is then of interest to measure 
the total number of events that contain back-to-back jets. 
Let Rjj/j(pT) be the ratio of events with back-to-back 
jets to the total number of events containing any jets, 
where Pt is the minimum value of px that any parti- 
cle must have to be counted as part of a jet. Thus, for 
example, for p T = 5 GeV, Rjj/j(pT) refers to the frac- 
tion of events containing two particles with pt > 5 GeV 
in opposite directions out of all events with any particle 
having pr > 5 GeV. Rjj/j (Pt) is then a measure of the 
fraction of space near the surface that can give rise to 
particles at high px through thermal-shower recombina- 
tion. It can have dependence on centrality. The precise 
relationship between £ and Rjj/j is a separate problem 
worthy of detailed investigation, especially if the experi- 
mental determination of Rjj/j is forthcoming. 



In a Monte Carlo calculation, such as that employed in 
, where space-time trajectories of the hard partons are 
tracked, it is possible to compute the parton momenta af- 
ter energy losses are taken into account. Those emergent 
partons then generate shower partons whose recombina- 
tion with thermal partons presumably results in a yield 
that can check the value of our mean suppression fac- 
tor £. In our treatment in the momentum space only, 
we determine £ by fitting the pion data, but once fixed 
the relative magnitude between TS and SS is also fixed. 
Moreover, there is no more freedom in the determination 
of the inclusive distribution of other particles, such as 
kaon and proton. 

D. Kaon production 



For the production of kaon the RF has the explicit 
form that follows from Eq. J5J with a — 1 and b = 2 |l3| 

2 3 

R K (p 1 ,p 2 ,p) = 12^ S(p 1 +p 2 -p) • (28) 
P 

Using this in Eq. JIJ yields for K + 



dN K+ 
palp 
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rP 

/ dp 1 p 1 (p-p 1 ) F U s(Pi,p-pi) • (29) 
Jo 



Note that because the kaon is not as tightly bound as 
pion, the RF is not as broad as that for pion, with the 
consequence that the factor pi(p — pi) 2 in the integrand 
forces the u and s quarks to have closer momenta than 
those in Eq. @J. 

There are four terms for F uS as in Eq. ©. The cal- 
culational procedure is basically the same as for pion. A 
slight complication arises from the strange quark being 
different from the light quarks. In the thermal compo- 
nent we shall simply attach a multiplicative factor A s for 
the s sector 



% = X s T , 



(30) 



where we set A s to be the Wroblewski factor at \ s = 0.5 
|19|. The thermal-shower recombination now has two 
terms 



J 



T( Pl )S s (p 2 )+S(pi)T s (p2)=ZC 



dkkfc(k) 



pie 



-Pi IT S S (P2 

k 



(31) 



The (SS) i and (SS) 2 contributions are as expressed in 
Eqs. (fT?)| and fHfy . respectively, with {j,f} identified as 
{u,s}. 

The four contributions to the kaon spectrum are shown 
in Fig. 4. No adjustable parameter has been used be- 



yond what has already been determined in the previ- 
ous section. The agreement between the sum (solid 
line) and the data is evidently very good. The data are 
for K® at 0-5% centrality and extend to px ~ 6 GeV 
poj . farther than for K+ [lElj. As in the case of pi- 
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g B = [B(a + l,/3 + j + 2)B(/3 + 1,7 + 1)]' 



(35) 



4 6 

P T (GeV) 



FIG. 4: Transverse momentum distribution of K + in Au- 
Au collisions. Data for K° are preliminary from |20B . The 
symbols for the four contributions are the same as in Fig. 1. 



ons the thermal-shower recombination is more important 
than both thermal-thermal and shower-shower (1-jet) for 
3 < pt < 8 GeV. As in Fig. 1 the shower-shower (2-jet) 
curve is shown for 8 y< jj = 0.01 and is negligible. There 
are some small differences in the various components of 
the pion and kaon spectra, but on the whole the two are 
basically similar. 



III. BARYON PRODUCTION BY 
RECOMBINATION 

A. General considerations 

The extension of the treatment in the preceding sec- 
tion to baryon production is conceptually straightfor- 
ward. The generalization of Eq. QJ is clearly 



dN E 



P- 



dp 



dpi dp 2 dp 3 
Pi P2 P3 



F(Pl,P2,P3) Rb(pi,P2,P3,p) (32) 



where F{pi,p 2 ,ps) is the joint distribution of three rele- 
vant quarks to form the baryon B. The RF is related to 
the non-invariant valon distribution by 

Rb(pi,P2,P3,p) = 5st 2/12/22/3 G B (yi, 2/2, 2/3), Vi = Pi/&3) 

where g s t is a statistical factor, and G B (yi, 2/2, 2/3) nas 
the general form 



For proton, y\ and y 2 refer to the momentum fractions 
of the U valons, and y 3 to that of the D valon. The 
exponents a, /3, and 7 have been determined in [22| to be 

a = £=1.75, 7 = 1.05 (36) 

from parton distribution functions of the proton. Note 
that although the U and D valons have the same con- 
stituent quark masses, they have different average mo- 
mentum fractions calculable from Eq. I|34|) 22] 

{y) v = 0.3644, (^ = 0.2712, (37) 

which satisfies the sum rule 



2 (y) u + (y) 1 



1 



(38) 



For hyperons, the lack of information about their par- 
ton distribution functions deprives us of any such detail 
knowledge of their valon distributions, and consequently 
of their RF's. Nevertheless, on the basis that the con- 
stituent quarks in all baryons are not tightly bound, we 
expect the exponents a, (3 and 7 for hyperons to be also 
in the range between 1 and 2. 

The 3-quark distribution now has more terms in the 
various possible contributions from the thermal and 
shower partons. Schematically, it takes the form 



qq'q" 



= TTT + TTS + T(SS) X 
+T(SS) 2 + (S(SS) 1 ) 2 + 



b {SSS) 1 
(SSS) 3 . 



(39) 



G B {yi, 2/2, 2/3) =9b y?y 2 3 yjS( yi + y 2 +y 3 -l) , (34) m-^B 



They are arranged in increasing order of the number of 
hard partons involved: the first term has none, the next 
three one, the following two two, and the last term three. 
We shall consider only the first four terms, since they 
involve thermal partons and the shower of only one hard 
parton. The fifth term involves partons from two over- 
lapping jets, and is ignored at RHIC energy despite the 
enhancement by T. 



B. Proton production 

Let us focus our attention on proton production at high 
Pt- As before, we omit the subscript T in referring to 
momenta in the transverse plane, and obtain from Eqs. 



dN v 



9 P 9st 



pdp p2a+ 7 +4 



dpi 



dp 2 (jPxPnYip-pi -P2) 1 F(pi,p 2 ,p-pi -p 2 ) 



(40) 
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The TTT contribution can be computed analytically. sideration is taken into account [lj . Equation Q41[l is not 

Using Eq. Q for each T, we obtain for the thermal spec- reliable at very small p, since the proton mass effect in- 

trum validates our essentially scale-invariant formulation for 

relativistic particles. Furthermore, at small px the prob- 

dN P = ° 3 pc -p/T B ( a + 2 > 7 + 2).B(a + 2, q + 7 + 4) i em becomes 3D, so our treatment should be modified 

pdp 6 B(a + 1,7 + l)S(a + 1, a + 7 + 2) ' accordingly. 

where C and T are given in Eq. i|10fl . and a, 7 in (|36|1 . 

The statistical factor <7 st is 1/6 when the spin- flavor con- The TTS contribution is 



dN TTS a C 2 F r f 

—^p- = epL-py+4 J dp 1 dp 2 (p lP2 r +1 (p-p 1 -p 2 y +1 Y / J dkkf t {k)u t {k lPu p 2 ,p) , 



(42) 



where 



_L p -(P~Pl)/T qu fPl\ P -(P-P2)/T Ou(P2 



Pi 

The TSS contribution is 
dN^ ss 



k ) p 2 



I ', ( k . p \ . p 2 . p) : — , "" ',s;(ii;..--e-WS*'(|)+ ; , ' .s;' ( J 1 ' ) 



-(pi+P2)/t od ( P-P1-P2 



P-P1-P2 



pdp 6p 



J d Pi d P2 (PiP2) a (p - Pi - P2) 1 J dkk f l (k)V i {k,p 1 ,p 2 ,p) , 



(44) 



where 



P - Pi - Pa 



^(,,Pi,P 2) p) = Pie-/, , f (m,)} +p2e -^ ^(i ^ ^ 

+ (p-pi - P2 )e-^-^)/ T {s« (|) , 



(45) 



r 



The curly brackets in Eq. I|45|l denote symmetrization as 
in Eq. (|17|) . Finally, we also have SSS which is sim- 
ply related to the fragmentation of a hard parton into 
proton. The latter has already been studied in 0- The 
corresponding FF, D\ (z), can be used here as in Eq. ijT§|) 
to give 



pdp p ^—^ 



dkmD^Q . (46) 



The result of our calculation for the four types of con- 
tributions are shown separately in Fig. 5. Their sum is 
shown as solid line, and agrees well with the data from 
PHENIX in Au-Au collisions at ^JJWn = 200 GeV for 
0-5% centrality We emphasize that there are no free 
parameters to adjust to achieve the good agreement. The 
result below px = 2 GeV is not compared with data, since 
the proton mass effect becomes important there and in- 
validates our scale-invariant formulation. All three con- 
tributions that involve thermal partons have about the 
same strength at pt ~ 4 GeV. What dominates at higher 
Pt is the TSS component until pt > 9 GeV where the 
SSS component takes over. The data available now are 
only for px up to 4.2 GeV and therefore cannot check our 



prediction in detail. At this point we can only conclude 
that the departure from the exponential behavior in the 
data for pt > 3.5 GeV is well accounted for by the two 
types of thermal-shower recombination, but not by di- 
rect hard-parton fragmentation. That may be regarded 
as empirical support for the role of thermal-shower re- 
combination. 

Having obtained both the proton spectrum here and 
the pion spectrum in Sec. 2, we can now calculate the 
p/ir ratio and compare it to the data. Since the calculated 
result shown in Fig. 5 is only for px > 2 GeV/c due to 
our neglect of the proton mass m p in our formulation, 
the p/n ratio that is represented by the solid line in Fig. 
6 is not shown for px < 2 GeV/c. It has the broad 
feature that the ratio is greater than 1 at px ~ 3 GeV/c. 
For low pt where the proton mass effect is important we 
adopt the Ansatz by replacing the factor 

which arises from Gb in Eq. J33J|, by m T ^ 2a+1 \ 
where tot = {rn^+p 2 ) 1 / 2 . The result is shown as dashed 
line in Fig. 6 and exhibits excellent agreement with the 
data 0, |23|. It should be noted that the amendment 
is at best kinematical, since no dynamical effect at low 
Pt has been considered. What is noteworthy is that the 
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p*(GeV) 8 



FIG. 5: Transverse momentum distribution of proton in Au- 
Au collisions. Data are from [lr|. The solid line is the sum 
of four contributions to the recombination of partons: TTT 
(dashed line); TTS (line with crosses); TSS, one thermal 
parton with two shower partons in one jet (dashed-dot line); 
SSS, three shower partons from one jet (line with squares) . 



maximum exceeds 1 at the peak, which is the anomaly 
that the fragmentation model cannot explain. The slow 
decrease of the ratio as px increases is a prediction of 
our model that can be checked by future data at px > 3 
GeV/c. 



1.5 



• PHENIX 

— proton mass neglected 

- - m T used 




FIG. 6: Comparison of calculated p/n ratio with data from 
Q. The solid line is the ratio of the solid lines in Figs. 5 and 1, 
without the mass effect taken into account. The dashed line 
is the result after p is replaced by tut in the way described in 
the text. 



IV. SAME-SIDE CORRELATION 

We have stated that because of the thermal-shower re- 
combination the structure of jets in AA collisions must be 



different from that in pp collisions. One way to probe that 
difference is to study same-side correlation of particles. 
In (3 we have calculated the 2-pion correlated distribu- 
tion in a w-quark initiated jet. That type of calculation 
can be incorporated into the study of same-side correla- 
tion in pp collisions by integrating over hard-parton mo- 
mentum and summing over all types of hard partons. For 
AA collisions one would consider 2q and 2q from thermal 
and shower partons and then recombine them to form 
two pions. Quantitative computation of the process will 
not be attempted here. Instead, we make some qualita- 
tive remarks in light of some data from RHIC that has 
some relevance. 

In a quantity Iaa is defined to be the ratio of all 
charged particles (above a modeled background) within 
a range of azimuthal angle A</> around the trigger par- 
ticle in Au-Au collisions to the same quantity (without 
background) in pp collisions. If there is no medium effect 
and if jets in A A and pp collisions are the same, then 
Iaa should be 1. It is found in ^3] that if the trigger 
momentum is in the range 4 < p^ lg < 6 GeV, the value 
of Iaa in central collisions is consistent with 1, though 
being more like 1.1. However, for 3 < p T rig < 4 GeV, 
Iaa is significantly higher, roughly 1.5 for V part from 
150 to 350. The former case for higher p^ ls has been 
regarded as evidence that the nearside azimuthal peaks 
in A A and pp collisions are similar [l?! 0] > but no ex- 
planation is given for the discrepancy when p^ ls is lower. 
The azimuthal distribution of the nearside particles is 
not shown in [HEil for 3 < p^ rig < 4 GeV. It is also not 
clear for that case whether the particles associated with 
the trigger are integrated over the range 2 < px < p^' s 
only, as it is stated explicitly for the case of 4 < p^ lg < 6 
GeV Hi. 



In our view there is no reason why Iaa should be equal 
to 1 for the same-side particles. Indeed, we regard the 
value Iaa ~ 1-5 for 3 < p T lg < 4 GeV to be an evidence 
in support of the enhancement due to thermal-shower 
recombination. The question that we should address is 
why Iaa is lower at higher p T ls . 

Although our formalism does not facilitate the calcula- 
tion of the azimuthal distribution, we can regard the inte- 
gral over A(j> under the same-side peak that enters the de- 
termination of Iaa as being roughly equal to an appropri- 
ate integral of the two-particle distribution dN/dPidP 2 
over Pi , with Pi being assigned the role of the trigger mo- 
mentum. Despite our restriction to collinear momenta, 
the integral J dP 2 accounts for all associated particles, as 
does / dA(j> in the data analysis by STAR |X3| - In our 
approach we have a hard parton at k creating a shower in 
which two partons at p\ and p2 combine separately with 
two thermal partons at p[ and p' 2 to form two pions at 
Pi = pi + p[ and P2 = P2 + p'2 ■ The distribution is 
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PlP2 J^k = J ^^^^(^)i(pi,^)r(pi)T( P ^(p lj y i ,P 1 )ii 7r ( ? > 2l p^P 2 ) . (47) 



For pp collisions the thermal partons are replaced by 
other shower partons initiated by the same hard parton 
that produces the ones at p 1 and p 2 , i.e., in Eq. (|47|l we 
make the replacement 

(SSU Pl , P2 )T(p' 1 )T(p' 2 ) -> (SSSS) 1 (p 1 , P2 ,p , 1 ,p , 2 )/^8) 

where (SSSS) i is a generalization of (SS) i in Eq. i|ltj|) 
to 4 shower partons. 

We now see that in pp collisions when the trig- 
ger momentum, Pi, is high, it forces pi or p[ in 
(SSSS)i(pi,p2,p'i,p' 2 ) to be high, with the consequence 
that P2 and p' 2 must be low. Shower partons Sf (z) at 
small momentum fraction z are known to have high den- 
sity. They are shown in Fig. 2 of Ref. 0. Indeed, the 
2-pion distribution in a hard parton based on (SSSS)i 
has been calculated in , where it is shown that the dis- 
tribution for the momentum fraction X 2 = P 2 /k of the 
(non-trigger) second particle becomes very high at low 
X 2 . The high density at low p' 2 renders (SSSS) i to be 
of the same order as T(p 2 ) in AA collision, resulting in 
an increase of the irn distribution in pp collisions. For 
that reason Iaa is lower (near 1) at higher p^ lg . In other 
words, since Iaa is dominated by the large number of 
particles in the low px range in the jets, it masks the dif- 
ferences in the structures of the jets (at intermediate and 
high pt) produced in nuclear and hadronic collisions. 

The drawback in using Iaa as a measure of same-side 
correlation is that it involves integrations in both Acj) 
and P 2 , a procedure that is likely to suppress the dis- 
tinctive features of thermal-shower recombination. It is 
recommended that the momenta of all particles in a jet 
are projected along the trigger momentum, Pi, and then 
the distribution in those projected momenta, P 2 , is de- 
termined for several values of Pi . 

V. CONCLUSION 

By showing the importance of considering shower par- 
tons created by hard partons, we have called into ques- 
tion the conventional paradigm in particle production at 
high px in heavy-ion collisions. The usual approach is to 
regard such particles as the products of parton fragmen- 
tation. We have shown that all particles are the result of 
parton recombination, including but not limited to the 
ones usually regarded as fragments. The important in- 
put that makes feasible our approach based entirely on 
recombination is the shower parton distributions derived 
in Ref. 0. Those distributions are determined by an- 
alyzing the fragmentation functions for parton-initiated 
jets in the recombination model. Once they are known, 



I 

it is conceptually unavoidable to consider the recombi- 
nation of thermal and shower partons in heavy-ion col- 
lisions. The result of such subprocesses turns out to be 
very important in the 3 < pt < 8 GeV range, as we 
have shown. The usual subprocesses of hard scattering 
followed by fragmentation are found to be unimportant 
until higher px- That results in a paradigm shift that 
has far reaching consequences. 

The phenomenon of energy loss of partons traversing 
dense medium can be related to experimental observables 
only by means of some valid model of hadronization that 
connects the partons to hadrons. If fragmentation is not 
important in the region of px under investigation, then 
serious modification of the quantitative implications of 
jet quenching must be considered, since the content of a 
jet has been altered from that produced in pp collisions. 
What such a modification should be has not been studied 
in this paper. We have only used an effective parameter 
£ to account for the fact that not all hard partons can get 
out of the dense medium to hadronize. We admit that 
such a method of treating energy loss is very rudimen- 
tary in this first attempt to study the effects of shower 
partons. Yet we have found consistency in being able to 
reproduce the spectra of pion, kaon, and proton. A more 
detailed investigation that tracks the space-time history 
of the produced partons would require an evolution code 
that is beyond the scope of this paper. 

Particles that are formed by thermal-shower recom- 
bination are part of a jet produced in heavy-ion colli- 
sions but are not present in jets in pp collisions. Thus 
the structures of jets produced in A A and pp collisions 
are different. To find evidences for such differences is of 
paramount importance in both experimental and theo- 
retical work to follow. The ratio Iaa discussed in Sec. 4 
on same-side correlation reveals a limited glimpse of that 
difference. We have suggested that the measurement of 
2-particle inclusive distribution at high px would provide 
a more accurate description of the jet structure and can 
be used to check the prediction that can be made in our 
framework. 

At this point the single-particle spectra that we have 
considered provide sufficient encouragement from the 
agreement with existing data to suggest that the recom- 
bination approach has captured the essence of hadroniza- 
tion at any px and that shower partons play an important 
role in the process. More data at higher px and on other 
species will give more stringent tests that the recombi- 
nation model must pass in order to establish the solidity 
required for a reliable mechanism of hadronization. 
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